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Chapter 5 ® 
Deriving the Lagrangian Density aaa 
of an Electromagnetic Field 


As part of understanding relativistic particles and electromagnetic fields, we need to 
have some awareness of the kinematic special theory of relativity, and then turn our 
attention to the dynamic aspect of charged particle motion in external electromag- 
netic fields. Thus, in this chapter we take a Lagrangian approach to the equations of 
motion and deal with the total energy involved with the motion of a particle. We also 
introduce the Hamiltonian in relation to the total energy of particle motion; analo- 
gous to classical mechanics, it relates to the corresponding kinetic and potential 
energy of the particle or system of concern. The Lagrangian approach is based in 
electrodynamics; thus, equations of motion are presented mainly as an avenue to 
introduce the concept of a Lorentz invariant action to a Hamiltonian, with the 
definition of the canonical momentum discussed in this chapter. 


5.1 Introduction 


Lagrangian field theory is a formalism in classical field theory. It is the field- 
theoretic analogue of Lagrangian mechanics, which is used for discrete particles 
each with a finite number of degrees of freedom. Lagrangian field theory applies to 
continua and fields that have an infinite number of degrees of freedom. 

The relativistic aspects of Hamiltonian and Lagrangian theories, which are rather 
abstract constructions in classical mechanics, have a very simple interpretation in 
relativistic quantum mechanics. They both are proportional to the number of phase 
changes per unit of time. The Hamiltonian runs over the time axis, whereas the 
Lagrangian runs over the trajectory of the moving particle, the t -axis, as illustrated in 
Fig. 5.1. 

The figure shows the relativistic de Broglie wave in a Minkowski diagram. The 
triangle represents the relation between the Lagrangian (£) and the Hamiltonian 
(J€), which holds in both relativistic and non-relativistic physics; it can be written 
with classical mechanics notation: 
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Fig. 5.1 The Hamiltonian and Lagrangian depictions scenario 


L=pv-H (5.1) 


The Hamiltonian counts the phase changes per unit of time on the vertical axis, 
whereas the term pv counts the phase changes per unit on the horizontal axis; v is the 
distance traveled per unit of time, whereas p is proportional with the phase changes 
per unit of distance—therefore the term pv. We now can understand the classical 
relation with g = x = v and write the following relation: 


a 


Og = (5.2) 


Thus, for the free classical relativistic particles we have the Hamiltonian energy and 
the po as follows: 








and pv = (5.3) 


Calculating the Lagrangian, we can see that the Hamiltonian is proportional to y, 
whereas the Lagrangian is proportional to 1/y, 





ao 2 
L = —(H — pv) = cama ee ( 1 )me (5.4) 


This is what we expect from time dilation. The moving particle has fewer clock-ticks 
by a factor y because of time dilation; we now check that: 


o£ oO v\ sl mm | 
ara ( 1 z)me\ = =a (5.5) 
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We have not yet discussed potential energy, which we need to expand on. To 
obtain the equation for motion of a relativistic particle in a potential field, we have to 
add the potential energy term V(q). In the non-relativistic case we have 


J =T(q) + V(q) 


5.6 
£=T(4) - Via) me 


1 
where T(q) = 5mnu is kinetic energy. 


The relativistic Hamiltonian and Lagrangian we have discussed, however, also 
include rest mass energy. The rest mass energy can be considered part of potential 
energy. Kinetic part T in the relativistic case can be obtained as follows: 

-1 


a oe re i 
H+L=2T=pv>T zpv = 5m 1-— ~ 5m (5.7) 








The approximation in this equation is valid when v < c. Using the term£ = 5m v7 in 
Eq. 5.5 provides the result for p = mv in the non-relativistic momentum, the same as 


we got here. 


5.2 How the Fields Transform 


As part of the derivation of the Lagrangian density of electromagnetic (EM) fields, 
we need to have some fundamental understanding of what field transformations and 
field tensors are all about; these two subjects are presented in this and the following 
sections. What we need to get our mind around is the general transformation rules 
for EM fields so that we can answer, for example, this question: 


1. Given the fields in system S, what are the magnetic fields in system Ss? 


To answer this question, our first guess might be that electric field Eis the spatial 
part of one 4-vector, and the magnetic field B is the spatial part of another. If this is 
true, then our intuition of one person’s electric field as another’s magnetic field is an 
incorrect assumption. Thus, let us begin by making an explicit assumption that we 
very implicitly can express in relativistic electrodynamics—that is, charge itself is an 
invariant case. 





Covariance | The property of a function to retain its form when the variables are linearly 
transformed ....under a covariance condition. 





Invariance | In mathematics an invariance is a property held by a class of mathematical objects 
that remains unchanged when transformations of a certain type are applied to the 
objects. 





Invariant A function, quantity, or property that remains unchanged when a specified trans- 
formation is applied. 
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+O Vv 


(c) 


Fig. 5.2 Illustration of moving capacitor 


This includes like mass, but unlike energy; the charge of the particle is a fixed 
number, independent of how fast the particle happens to travel from one point to 
another. Thus, we assume that the transformation rules are the same no matter how 
the fields are produced—that is, electric fields get generated by moving currents and 
transform the same way as those set up by stationary charges. If this is not the case, 
then we would need to abandon the field formulation altogether because what it is 
essential for a field theory is that the fields at a given point provide all the 
information one needs to know. Yet, in case of EM fields, one does not have to 
have extra information, electromagnetically, regarding their source [1]. 

Bearing the preceding argument in mind, we can consider one of the simplest 
possible electric fields, which is uniform in the region between the two-large parallel- 
capacitor plates, as shown in Fig. 5.2a, by assuming that the capacitor is at rest in So 
and carries surface charges +o9. Then, we can write the electric field Eo in terms of 
surface charge density oo in free space and €) = 8.85 x 107"? C?/N x m’, this being 
the permittivity of free space, as well as a unit vector ¥ in the y-direction, as follows: 





pee (5.8) 


€0 
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If we consider the same capacitor from system S moving to the right at speed vp, as 
shown in Fig. 5.2b, however, we can see that these system plates are moving to the 
left, but the field takes the form of the following relationship [2]: 


E=—3 (5.9) 


The only difference between Eqs. 5.8 and 5.9 is the value of the surface charge o, 
but we need to see whether that is the only difference. Equation 5.9 is written for a 
parallel plate capacitor, which came from Gauss’s Law, whereas Gauss’s Law is 
perfectly valid for moving charges, where this application totally relies on the 
symmetry issue of the problem. If the plates tilt, however, as we can see in 
Fig. 5.2c, is it possible to say that electric filed E still is perpendicular to the moving 
direction of the plate? The answer is that it does not matter, and the field between the 
plates is the creation of the superposition of the surface charge of +o. The one from 
the —o nevertheless would run perpendicular to the plates that are in motion, even in 
a tilted condition. For the surface charge, —o, field would orient as indicated in 
Fig. 5.2c; changing the sign of the charges reverses the direction of the field, and the 
vector sum kills off the parallel components [2]. 

Now we are at the point where we can consider the total charge on each plate and 
claim that as invariant and the width (w) is unchanged; however, the length (1) is 
Lorentz contracted by a factor of: 


1 
—=,/1-v2/e 5.10 
rma o/ (5.10) 


So that the charge per unit area is increased by a factor of 7p as indicated here: 
o = 7000 (5.11) 


Thus, accordingly, we can say that: 
B= Yokg (5.12) 


The symbol is clearly an indication of a perpendicular component of electric 
field E, while symbol || is a presentation of a parallel component of same field in the 
direction of the motion of system S, considering that the capacitor lined up with the 
yz-plane, as shown in Fig. 5.3. At this time, it is plate separation d that is Lorentz 
connected, whereas / and w, consequently, o, are the same in both frames. Because 
the field does not depend on d, it follows that: 


~ =) 
What we have learned so far from the previous discussion we can put into 


perspective with the following examples given by Griffiths [2]. Here, we are going 
to show a simple problem of an electric field with a point charge in uniform motion. 
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Fig. 5.3. Two parallel y 
charged plates 





Example 5.1 A point charge gq is at the origin in system So. 


Question What is the electric field of this same charge in system S, which moves to 
the right at speed vp relative to system So? 


Solution In system So, the field is: 
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Are ro 
1 GZ 
Ea Ameo (+2 2. ,2)3/2 
0(x r Yo %) 
1 4X0 
Exo 3 
; 3/2 
4ne0(x3 +33 + 8) 
1 qY0 
En = 
© ARoGg +95 +4)" 











From the transformation rules of Eqs. 5.12 and 5.13, we can state that: 














ee 1 qXo 
<a 4m€o (42 4 y2 1. 2)3/2 

(3 +93 +2) 

YoIo 
EB) = YoE\o = 3/2 
Ane (xs +3 +28)" 

1 709% 
E. = YoE 20 = 3/2 
Anco (x5 +9 +)" 
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Fig. 5.4 Illustration of 
point charge in system So 








The expression is still in terms of the system Sp coordinates (Xo, yo, Zo) of the field 
point (P), and preferably we need to write them in terms of the system S coordinate 
of point (P). From the Lorentz transformations, or actually, the inverse transforma- 
tions point of view, we can write: 

Xo = Yo(x + vot) = 7oRx 

20: y =R, 

Zo = Zz =R, 


where R is the vector from charge g to point P, as depicted in Fig. 5.4. Thus, we can 
write: 





aa 1 vyogR 
Ane | 2 R2 2 2 in 2 3/2 
7jR? cos?(0) + R? sin ?(0)] 
pa (5.14) 
I q(i-v/e?) oR 





~ Ane (1 — (v3 /c?) sin2(0)]°” R? 


This equation then is identified as the field charge in uniform motion; thus, we get the 
same result if we use a retarded potential equation, known as Lienard-Wiechert 
potential for point charges if they move on a specific trajectory [2]. 

This is a significant and efficient derivation, and shows that the field points away 
from the instantaneous as opposed to the retarded position of the charge; E,, gets a 
factor of yo from the Lorentz transformation of the coordinates, and E, as well as E, 
pick up their factors from the transformation of the field. It is the balancing of these 
two yos that leaves E parallel to R. 

Nonetheless, Eqs. 5.12 and 5.13 are not the most general transformation laws that 
we can express because we start with system Sp in which the charges were at rest and 
where, consequently, there was no magnetic field. For us to derive the general rule, 
what we need to start with is out in a system with both electric and magnetic fields 
simultaneously. For this purpose, S itself will serve very nicely. 
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In addition to the electric filed, 
Ey = (5.15) 


there exists a magnetic field because of the surface currents, as shown earlier in 
Fig. 5.2b: 


> 


Kz = Fouot (5.16) 











Using the right-hand side (RHS) rule, as it was mentioned in Chap. | of this book, 
this field points in the negative z-direction and its magnitude is given by Ampére’s 
Law as: 


B, = —[lyovo (5.17) 


Using a figure, such as depicted Fig. 5.5, for a third system S, traveling to the right 
with speed v relative to S, the field would be: 


Ey = B, = —[pjov (5.18) 


where 0 is the velocity of system S relative to system So, and it can be written in the 
following form: 


= v + 00 
v= 
1 + vv9/c? 
1 (5.19) 
y= /1 — 0? /c? 
and 
6 = Yoo (5.20) 


Fig. 5.5 Illustration of = 
systems So, So, So and their ya 
motions v4 


You 


Vo v (v relative to S) 
=> => 
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The only thing that remains is to express electric filed E and magnetic field B, as 
shown in Eq. 5.18, in terms of F' and B (see Eqs. 5.15 and 5.17). In view of Eqs. 5.11 
and 5.20, we have the following: 


Be (z)= 
Yo/ €0 
B, = — (Z) Hood 
YO 


With a little manipulation of algebra, we find the following result: 


y Vive eae ~) (5.22) 


(5.21) 











Yo Vl-P/e2 /1—v?/e? f em 
where 
: 5,23 
t= fi-ve oi 
Thus, we can write: 
B= (1 Es = = (6, = Bs) (5.24) 
c’ / 9 C7 E0Ug 
whereas 
= vU0 Vv + V9 
B= (1 ai ) 
eee ee moe(, + roula) (5.25) 


= 7(B, — UeovEy) 
Or, because pio€o = I/c*, then Eqs. 5.24 and 5.25 reduce to the following form: 
Ey = 1 (Ey _ vB-) 
= ”) (5.26) 
B, = 7(B. > “Ey) 
a 


This equation indicates how E,, and B, transform—for E, and B, simply align the 
same capacitor parallel to the xy-plane instead of xz-plane, as illustrated in Fig. 5.6. 
Thus, the fields in system S are written: 


£0 (5.27) 
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Fig. 5.6 Two capacitors y 
parallel to the xy-plane 





Using the RHS rule of Chap. | again, we get B,. The rest of the argument is 
identical—everywhere we had EF, before reads E., and everywhere that we had B, 
reads —B,,; thus, we can write: 


E.= r(E: oF vBy) 
7 D (5.28) 
B, = 1(By + SE.) 


As for the x-components, we already have seen, by orienting the capacitor parallel to 
the yz-plane, that: 
E,= FE, (5.29) 


Because in this case there is no accompanying magnetic field, we cannot 
conclude and deduce the transformation rule for B,. Nevertheless, another configu- 
ration, such as the one in Fig. 5.7, which is an illustration of a long solenoid aligned 
parallel to the x-axis and at rest in system S. The magnetic field within the coils is 
then [2]: 


By, = eon! (5.30) 


where n is the number of turns per unit length and / is the current in the coil. In 
system S the length contracts, so n increases: 


n=yn (5.31) 


On the other hand, time dilates. Thus, the system S clock that rides along with the 
solenoid runs slow, so the current (i.e., change per unit time) in S is given by: 
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Fig. 5.7 Depiction of long y 
solenoid aligned parallel to 
the x-axis 
CUCLCCCCCUAN 7 
z 1 
- | 
T=-I (5.32) 


The two factors of y exactly cancel, and we conclude that: 
B, = By (5.33) 


Like E, component B, which is parallel to the motion, is unchanged [1]. 
Let us now collect all the components of transformation in order to have a 
complete set of transformation rules: 














E,=E, Ey=y(Ey—0B,)  E,=y(E,+0By) 
a _ v = Dv (5.34) 
B.=B, By= (By - <E,) B, = 7(B. = <E;) 
uae 2 
The following are two special cases that we need to pay close attention to: 
1. If B=Oin system S, then we can write: 
2 D px ~ Dive. hs 
B= Ya (Ey — Eyz) = a (E.¥ — E,2) (5.35) 
or, since V = vx, then 
> ly, .2 
Bs (¥ x E) (5.36) 
2. If E=Oin system S, then 
E = —yo(BS — By2) = —v(B5 — Bz) (5.37) 


or 


E=%xB (5.38) 
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5.3. The Field Tensor 


As part of electrodynamics and relativity, we can expand our knowledge into the 
field tensor, where the classical electrodynamics already is consistent with special 
relativity that can be extended and be applied to Maxwell’s equations and the 
Lorentz force law It is possible, legitimately, in any inertial system, even in an 
invariant Lorentz transformation, to show how the field transforms. See Appendix B 
for more details [2]. 

To understand what we discussed as part of developing Eq. 5.6 and the scalar 


(I E I" _ *||B\)’). one of the two Lorentz invariants, we consider the following 


tensor field analyses; refer to Appendix B for more essential information and 
derivations. Still, bear in mind that a 4-vector transforms by the rule shown in the 
following equation: 


a” = AH’ (5.39) 


Note that the summation over v implies that a A is the Lorentz transformation 
matrix. If S is moving in the x-direction at speed V, then A has the form, 


y vB 0 0 
= 0 0 
a2 7 (5.40) 
oo 9 0 0 
0 0 00 


A‘ is the entry in row y and column v. A second-rank tensor is an object with two 
indices that transforms with two factors of A, which is one for each index as 
demonstrated in the following equation, where we have: 


eh (5.41) 


A tensor in four dimensions has 4 x 4 = 16 components, which we can display in 
a4 x 4 matrix array: 


f00 {01,0203 
10 ll 41213 
t t t t 
wv 
a 20 21 722 423 (5.42) 


po p 1 p2 3 


But, the 16 elements need not all to be different. For instance, a symmetric tensor has 
the property as follows: 


t” =f" (Symmetric Tensor) (5.43) 
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In this case, there are 10 distinct components and 6 out of 16 represent the repeats 


(91 = 710, 0? = 70 3 = PO fl? = P13 = A! and 73 = £). Similarly, and 


antisymmetric tensor obeys the following relation as well: 
#” = —f'¥ (Antisymmetric Tensor) (5.44) 


Such an object has just 6 distinct elements of the original 16 components, where 
6 are repeats of the same ones as before, only this time with a minus sign in front of 
each element, and 4 of them are zeros; they are ic f, r?, and ry), Thus, the 
general form of an antisymmetric tensor has the following format: 


0 pl 792 13 


, 0 4 712 413 
= 20 21 23 (5.45) 
—t —t 0 " 


po _f1 _22 9 


Now let us see how the transformation of the rule defined by Eq. 5.41 works for 
the 6 distinct components of an anti-symmetric tensor. Thus, starting with compo- 
nent —t'°, we have the following relation: 


Sana (5.46) 


According to Eq. 5.40, a = 0; however, unless 4 = 0 or 1, and Sa = Ounless o = 0 
or 1. Thus, there are four terms in the sum: 


—1°! = ApAgt™ + AQAIO! + AVAGt? + ALAGT (5.47) 


On the other hand, °° = r'! = 0, while pi = a then we can write the 
following relations: 


=11 = (AGAY — APAG) = [7 - (7p)] 0" = 2 (5.48) 
Working out the other components, similarly, we obtain the following results: 


=! = pl 902 = (12 — ppl?) 13 = (1 — Bp!) 
—3 = —p! = (8! — pr) 12 = y(r!2 — py®) 
These are precisely the rules we can derive on physical grounds for the EM field in 


Sect. 5.2 [1]. 
In fact, we can construct field tensor F“” by direct comparison as: 


(5.49) 


Ey E, E 
PF! = a F°2 =. PF = 

c ; c , © (5.50) 
Pea, Fak, Fae, 
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Written in Eq. 5.50 as a matrix array, we can produce the following: 


0 E,/c E,y/c E,/c 


-E,/c 0 B, —B, 

FY = (5.51) 
=E,J/c =—B, O B; 
7 nf C By —B, 


With the results obtained in Eq. 5.51, relativity completes and makes a perfect job 
by Oersted, combining the electric and magnetic fields into a single entity as F””. 
Note that some authors prefer the other notation convention—F?! = E,, F i cB,, 
and so on—and some use the opposite signs. Accordingly, most of the equations 
from here on will have a little different appearance, depending on the text you are 
reading. 

If we follow the argument that yields Eq. 5.51 results, with some cautious care, 
we can notice that there is a different way of relating an anti-symmetric tensor with 
electric field E and magnetic field B. Instead of comparing the first line of Eq. 5.51 
with the first line of Eq. 5.49, and the second with the second, we can relate the first 
line of Eq. 5.51 to the second line of Eq. 5.49 and vice versa. This leads to and results 
in the dual tensor, G“”, as follows: 


0 B BB, 
—B, 0 —E,/c —E,/c 
G! = : (5.52) 
—-B, E,/c 0 —E,/c 
=B,. -=Byje Hyfc 0 


G"” can be obtained directly from F“” by the substitution of E/c > Band B > E/c. 
Notice that this operation leaves Eq. 5.51 unchanged; this is why both tensors 
generate the correct transformation rules for E and B. 


5.4 The Electromagnetic Field Tensor 


The transformation of electric and magnetic fields under a Lorentz boost was 
established even before Einstein developed the theory of relativity. We know the 
E-fields can transform into B-fields and vice versa. For example, a point charge at 
rest gives an electric field. If we boost to a frame in which the charge is moving, there 
is an electric and a magnetic field. This means that the E-field cannot be a Lorentz 
vector. We need to put the electric and magnetic fields into one tensor object to 
properly handle a Lorentz transformation and to write the equations in a 
covariant way. 
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The simplest way, and the correct way, to do this is to make the electric and 
magnetic fields components of a rank 2 anti-symmetric tensor as follows: 


0 B. —By —iE, 


-B, 0 B, —iE, (5.53) 


The fields can be written simply in terms of the vector potential, which is Lorentz 
vector A, = (A, ip), thus Eq. 5.53 can be expressed in the following form: 


_ OA, OA, 
Saye: ee 








Fuy (5.54) 

Note that this is automatically anti-symmetric under the interchange of the 
indices. As before, the first two source-less Maxwell equations are automatically 
satisfied for fields derived from potential. We may write the other two Maxwell 
equations in terms of 4-vector j,, = ( J icp) in the following form: 








OF yw Ju 
ie (5.55) 
Comparing Eqs. 5.54 and 5.55, we can write the following: 
0 (dA, OA\ iy 
= 5.56 
Ox, & =) Cc ( ) 


Of course, we have not yet quantized the theory in Eq. 5.56. 

For some peace of mind, let us verify a few terms in the equations. Clearly, all the 
diagonal terms in the field tensor are zero by anti-symmetry. Taking some examples 
of off-diagonal terms in the field tensor, check the old definition of the fields in terms 
of potential: 


B=VxA 

» «, 10A 
B= 25 
Foe (Vx A) =8, 
Fy = SB AL (Fx A), =—2, 





_ 0A; OAg 10A; O(ip) —.(10Ai  OP\ (0g , 1 0A; 
Ox, Ox; ic Ot Ox; NC Ot * Ox;) : ons Ot 
(5.57) 
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Let us also check what the Maxwell equation says for the last row in the tensor, as 
follows: 











OF 4, _J4 

Ox, C¢ 

OF 4 | icp 

Ox; oc 

(iE) _ (5.58) 
Ox; 

aE 

Ox; Pp 

V-E=p 


We will not bother to check the Lorentz transformation of the field here. 


5.5 The Lagrangian and Hamiltonian for Electromagnetic 
Fields 


There are not many ways to make a scalar Lagrangian from the field tensor. We 
already, from Eq. 5.55, know that 


OF w _ Ju (5.59) 


Ox, Cc 





In addition, we need to make our Lagrangian out of the field, not just the current 
again; x, cannot appear explicitly because it violates symmetries. Also, we want a 
linear equation, so higher powers of the field should not occur. Term mA,A" is a 
mass term and would cause fields to fall off faster than 1/r (i.e., we later show how 
the term 1/r comes into play). Thus, the only reasonable choice is: 


F,)F“” = 2(B? — E’) (5.60) 
One might consider the following relationship: 


Ry A a (5.61) 


But that is a pseudo-scalar, not a scalar completely—that is, it changes its sign under 
a parity transformation. The EM interaction is known to conserve parity, so this is 
not a real option. As with the scalar field, we need to add an interaction with a source 
term. Of course, we know electromagnetism well, so finding the right Lagrangian is 
not really guess work. The source of the field is the vector ig thus, the simple scalar 
that we can write is j,A“. 
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The Lagrangian for classical electricity and magnetism we will try is: 
1 yl 
Lem = ——F yy FY + =I (5.62) 
4 c 


In working with this Lagrangian, we will treat each component of A as an indepen- 
dent field. 
The next step is to check what the Euler-Lagrange equation provides: 














0 o£ o£ 
ax (seax vag) ae Oe) 
where 
1 pw be gy 1 fA" BAK) (3A” BAM 
a= gee ae 4 (= =~) ( =) pe) 
o£ ee 0 OA’ 0A*\ (0A? OA* 
O(OA"/Oxy) 4. 0(0A“/Oxy)\ Ox, OXe] Ox, OX 








1 0 ( OA’ 0A° 5 OA’ a) 








40(0A"/0x,)\ Ox, Ox, Ox, OX5 (5.65) 
ae 0A’ OA" 
~ 4 \ Ox, Ox, 
= Fy, = Fy 
0 OL 
yy 
oa. A 
9 pw _Iu_g (5.66) 
Ox, c 
0 J 
yo _ Sn 
oa a 


Now that we have four independent components of A“ as autonomous fields, we 
have four equations or one 4-vector equation. The Euler-Lagrange equation gets us 
back to Maxwell’s equations with this choice of the Lagrangian. Clearly, this 
justifies the choice of £. 

It is important to emphasize that we have a Lagrangian-based, formal classical 
field theory for electricity and magnetism, which has the four components of the 
4-vector potential as the independent fields. We could not treat each component of 
F,y as independent because they obviously are correlated. We could have tried using 
the six independent components of the anti-symmetric tensor, but it would not have 
given the correct answer. Using the 4-vector potentials as the fields does give the 
correct answer. Electricity and magnetism is a theory of 4-vector field A”. 
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We also can calculate the free field Hamiltonian density—that is, the Hamiltonian 
density in regions with no source term. We use the standard definition of the 


Hamiltonian in terms of the Lagrangian: 








PPOs EL SN ae 
O(OA"/Odt)) Odt ~  \O(OA¥/Oxa)) Oxa 
We calculated before that: 


axe 
O(OAJox,) 


which we can use to get: 


o£ 


B(OA/Ox,) 1 


Thus, the free field Hamiltonian density becomes: 





Oe nae ee 
( 4) Ox4 
OAK 1 . 
a 
Therefore, 
OAK 1 
H = Fy4~—t+-F FF” 
Wf Ox4 ‘s 4 ee 











ae 


(5.67) 


(5.68) 


(5.69) 


(5.70) 


(5.71) 


We will use this once we have written the radiation field in a convenient form. In 
the meantime, we can check what this gives us in general in a region with no sources: 


OA* 1 
He= Fis Fag to—) t—F yy FY’ 
1 4 +) +3 m 
OA* 1 


= ya( Fa +z \+ ike” 
i Oxy 4 


Ae) 5 EY 
uw «4 ss 


= FP ay — Pua +o 








1 2 2 
-+5(8 E’) 








(5.72) 
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If we integrate the last term by parts (and the fields fall to zero at infinity), then that 
term contains a V - E, which is zero with no sources in the region. We therefore can 
drop it and are left with: 


d= ; (E* + B’) (5.73) 


This is the result we expected: the energy density and an EM field. (Remember the 
fields have been decreased by a factor of 4 compared to centimeter-gram-second, 
CGS, units.) 

We will study the interaction between electrons and the EM field with the Dirac 
equation later. Until then, the Hamiltonian used for non-relativistic quantum 
mechanics will be sufficient. We have derived the Lorentz force law from the 
Hamiltonian. 


5.6 Introduction to Lagrangian Density 


In considering the special relativistic EM field, we understand that assuming a 
Lagrangian density of the form results in: 


and following the Euler-Lagrange equation recovers Maxwell’s equation. In this 
equation F“” is the Maxwell field tensor and c is the speed of light, c* = l/ey (not 
necessarily in a vacuum). The four-current is re = (pe, j); the four-potential is A 
= (p /c, A); the Compton wave number for a photon with mass (m) takes the form 
k = 2amclh; and his Plank’s constant. 

In the following we will prove that a compatible Lagrangian density for the EM 
field in empty space is, 


ate) 1/2 
E\l” — c?||B es 


which will be proof of Eq. 5.74. In following section, we take a very general 
approach to the more classical electrodynamics of Lagrangian density for the EM 
field. 

Equation 5.75 is the Euler-Lagrange equation produced from this Lagrangian 
and are Maxwell’s equations for the EM field. Lagrangian density is derived by a 
trial-and-error procedure, not by guessing, and one also can find the Lagrangian in a 
more difficult and complicated case. (See the proof provided in Law [3], where the 
interaction between a moving mirror and radiation pressure is a Hamiltonian 
formulation.) 
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Maxwell’s differential equations of the EM field in empty space are: 


Vx B= - (5.76a) 
ee oe, 1 BE 
Vx 8=4)+-— 76b 
x B= Wi + 3a (5.76b) 
v= (5.76c) 
€0 
V-B=0 (5.76d) 


where E is the electric field intensity vector, Bis the magnetic-flux density vector, p 
is the electric charge density, and i is the electric current-density vector. All 
quantities are functions of the three space coordinates (x, x2,x3) = (x,y,z) and 
time t = x, for time coordinates of the fourth dimension. 

From Eq. 5.76d, magnetic-flux vector B may be expressed as the curl of vector 
potential A and is written as follows: 


B=VxA (5.77) 


In addition, from Eqs. 5.77 and 5.76a, we can conclude that: 
af =e. 0A 
V E+ |= 5.78 
« (2434) 5789 
So the parentheses term may be expressed as the gradient of a scalar function: 


~. OA 
Bag > -~V¢ (5.79a) 


Equation 5.78 can be rearranged as follows as well: 


> 


E= ae ae rs (5.79b) 


Therefore, the six scalar variables, the components of vectors E and B, can be 
expressed as functions of four scalar variables, the scalar potential ¢, and three 
components of vector potential A. Inserting the expressions of E and B, Egs. 5.77, 
5.79a, and 5.79b, respectively, into Eqs. 5.76b and 5.76c, we obtain the following 
result: 





a ee + 10 0A 
Vx (Ox A) =mi+a5( Vo =) (5.80) 
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and 





Wp (ae (5.81) 











Given that the vector identity is of this form: 
Vx(VxA)=V(V-A)-V7A (5.82) 


Per Eqs. 5.82, and 5.80, we find the following result: 





Oe 


0°A a _ 1d¢ : 
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5.7 Euler-Lagrange Equation of the Electromagnetic Field 


Now the main goal and task is to find a Lagrangian density £ function of the four 
“field coordinates” and their first-order derivatives as follows: 


L=L(api, Vj) j=1,2,3 (5.84) 


such that the four scalar EM fields, Eqs. 5.81 and 5.83, are derived from the 
Lagrange equations: 


k=3 
] Of | sn 8 |_ 0% oF 9 9493 (5.85) 


Ot 3 on; > Bx 3 On; On; 
Ot Oxy 
This can be simplified in notation to: 
Of 0£ = 
(32) oe 


Here the Lagrangian density £ is a function of the following: 








o£ 


0(Vnj) 


ax 
ah, pO 8 5.86 
On; J ( ) 











1. The four “field coordinates” 


ny = Aj (x1, x2, x3, t) (5.87a) 
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My = A2(x1,%2,%3, 1) (5.87b) 
13 = A3(x1,%2,%3, 1) (5.87c) 
Ng = Aa(x1,X2, 3, t) (5.87d) 

2. Their time derivatives 
= Mea, (5.88a) 
n= Jas es (5.88b) 
13= f =A a4, (5.88c) 
i a = = =¢ (5.88d) 

3. Their gradients 

Vn, = VAL, Vig = VA, Vin3=VA3, Vig = Vb (5.89) 


We express Eqs. 5.81 and 5.83 in forms that are similar to Lagrange Eq. 5.86: 





Feet, 
$(¥-A)+9-(96) -(-£) = (5.90) 











Finally, we can write: 


0 (OA , OP\ | (OA 

Or ( ar T x] T V- f (#- va‘ 
The Lagrange Eq. 5.86 for j = 4—that is, for 74 = @ is given as: 

0 (o£ = o£ o£ 

eee (peta V-|— =| -~ =0 5.92 
ai (ag) +? Lawes] ~ 34 ae 


Comparing Eqs. 5.90 and 5.2, we note that the first could be derived from the second 
if, 








4-0 (5.91) 














aL . o£ af  p 
SS yay yh. ee 
Og ave ee 


so that Lagrangian density £ must contain, respectively, the terms: 
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ae +s 1 
La, =(V- Ald, Lon = 5IIVAl?, La, = (5.94) 


Consequently, their sum becomes: 
- ~. 1 
Lu = La + Latte =(F Ab +zIvel?—-—% (5.98) 
0 


We suppose that an appropriate Lagrangian density £ would be of the form: 
L=Lat+L, (5.96) 


and because £, produces Eq. 5.90, we expect that £;—to be determined—will 
produce Eq. 5.91. This expectation would be correct if Eqs. 5.90 and 5.91 were 
decoupled; for example, if the first contains only ¢-terms and the second contains 
only A-terms. But, here this is not the case: £L, containing A-terms would participate 
to the production of Eq. 5.91; moreover, £, would participate in the production of 
Eq. 5.90. It is possible that this could mutually destroy the production of the 
equations as we had expected. Nonetheless, here we follow a trial-and-error proce- 
dure, which will direct us to the correct answer, as we can see in the following. 
Now Lagrangian Eq. 5.86 for j = k = 1, 2, 3—that is, for 4, = A,—is: 


0 ( o£ oL o£ 





Comparing Eqs. 5.91 and 5.97, we note that the first could be derived from the 
second if, 


dL . d¢ of (0A OL ik 
Fy Pen O(VA,) (72-va) ee (5.98) 


From the first part of Eq. 5.98 the £, of Lagrange density £ must contain the 
following terms: 


1). 2 O¢, 
5llAal| +a An k= 12.3 (5.99) 
So their sum with respect to k is: 
1) 34/2 % 
Lp, =5 A|| +Vq@-A (5.100) 














From the second part in Eq. 5.98, the £, of Lagrange density £ must contain the 
terms as follows: 


1,/0A 
2 OA A, — | VAgI? k=1,2,3 (5.101) 
2 | Oxe 
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So their sum with respect to k turns into this format: 





k=3 7 
La aie) a4 Td I| V Axl” (5.102) 
2 2 kel Ox, 








From the third part in Eq. 5.98, the £, of Lagrange density £ must contain the 
terms: 


JiAk 
€0 


k=1,2,3 (5.103) 


Therefore, their sum, with respect to k, becomes: 


4 (5.104) 


From Egs. 5.101, 5.103, and 5.104 the £, part of Lagrange density JL is: 
Lp =Lp, + Lp, + Lp, 


4 +v0 ee ie 0A 4 IV Aull? 
oa kel Oxy : : 


7 A 
ad (5.105) 
€0 





=; 
2 





Finally, from Eqs. 5.96 and 5.104 for densities Land L,, Lagrange density 
L=L, + Lp is given as: 





L=LatLy 
aes 1 
= (7 -A)b +51V0r -2 
2 £0 (5.106) 
Ly) > ||? 2 Ages OA a| oA 








Yet, this is a “wrong” Lagrange density; thus, we need to perform the following 
error-trial-final success procedure (see next section). 


5.7.1 Error-Trial-Final Success 


Insertion of this Lagrange density expression into the Lagrange equation with 
respect to @ (i.e., Eq. 5.92) does not yield Eq. 5.81 but gets: 


0 


2 
_V o-( 


2V - A) ete (wrong) (5.107) 
£0 
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The appearance of an extra (V . A) is because of the term (Vo -A ) of £, and that is 
why the Lagrange density given by Eq. 5.106 is not an appropriate one. To resolve 
this problem, we must look at Eq. 5.90 (i.e., Eq. 5.81) from a different point of view 
as follows: 


V-(Vp+4)- (-£) = (5.108) 


Comparing Eqs. 5.108 and 5.93, we note that the first could be derived from the 
second if in place of Eq. 5.93 we have the following result: 


2Voe4,. —e-—= (5.109) 





So in place of Eqs. 5.94 and Eq. 5.95, respectively, the new form of the equations is 
as follows: 


£,, = 0, £1, = IVA? + V6- A, Li =L.,=—— (5.110) 


1 S 
= 2, +01, + Li, = 5IIVOl? + (Vo A) — 2 (5.111) 


€0 
Now it is necessary to omit from L,, Eq. 5.100, the second term (Vo . A) because 
it appears in to (see the second term of Eq. 5.110). Thus, we have in place of 
Eq. 5.100, the following form: 
' 1 


n= 5 A (5.112) 














Although £,, and £L,, remain unchanged, as in Eqs. 5.102 and 5.104, we have: 














f.=2 aS oA VAx — ||VAal|? (5.113) 
By — P~Py = 2 A Ox, k k ‘ 
J-A 
3, = Lp, = . (5.114) 
In place of Eq. 5.105 we now can write: 
/ / / / 
Lg=L5, +L, + Lp, 
Z A (5.115) 





Ijsj2, Lows| OA , 
=+5||4| 422 VA, — || VAql2| +4 
2, me Oxx y all €0 
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Finally, for the new Lagrangian density we have the following in place of Eq. 5.106: 
L=L,4+2L, 


aig VieuaA = ee 
all ol + Vo- (5.116) 


Dl 

















Hey ‘fe -VAg = || VAul|? 





Density L of Eq. 5.116 is obtained from density £ of Eq. 5.106 if we omit the term 
(V- A). Thus, £ is independent of @. 

In the following equations the brace over the left three terms groups that part of 
density L', which essentially participates in the production of electromagnetic 
Eq. 5.81 from the Lagrange equation with respect to @, Eq. 5.92. Although the 
brace under the right four terms groups that part of density L', which essentially 
participates in the production of electromagnetic Eq. 5.83 from the Lagrange 
equations with respect to A,, Az, A3, Eq. 5.87. Thus, we have: 


with respect to 














: k=2 > 
slVal? - re Ve: A +54] eh wa (Vag? | = 
k=1 
= hive? 24 va-a+ dal +4 oe vay ae? | +L 
with respect to A 
(5.117) 


Note the common term (Vo . A). 
Recording the terms in Eq. 5.116 of density SL’, we have: 











‘|/a|| +2 elo gi I — 
= 54] + Ie? + V6-A—S2 SI VaeP — 5 Vac! +(-0b +7 A) 
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(5.118) 
that is, 
1 aA] 1 1 
a V 2 V A 2 Ra eae rae cf 
Lf =5|-Ve— Z| gel xP + (ob 47-A) 6.119) 
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or 


= ))2 = 1)2 
Ell — c?||B 
A -eair t 


x ; a pp+j-A) (5.120) 





Now if density £° must have dimensions of energy per unit volume, we define 
Lem = Eo ; therefore, we can write the following result: 























a2 9) By2 - 
Lem = & ie = Ml | po+j-A (5.121) 
keeping in mind that 
dal’ y=) ; 
ra" - |--+-F =A +iver+20ve-A) (6.1228) 
31/2 = =||2 = 2 OA 
aI? =F x Al => [ivan -$4 vay (5.122) 


The scalar (EI = |||") is one of the two Lorentz invariants (e.g., see 


subsection in the preceding description of the field); the other one is E - B. The scaler 
essentially is equal to constant time ¢,,e”", where e”” is the anti-symmetric field 
tensor (see next subsection as well). (More details can be found in Appendix A under 
the topic Relativity and Electromagnetism.) 

On the other hand, the scalar (—pp + j: A) is essentially the inner product J,,A“ 
of two 4-vectors in Minkowski space. The four-current density J4 = (cp, ri ) and the 
four-potential A” = (¢ Ic, A) is a Lorentz invariant scalar too. Therefore, Lagrangian 
density L,,,, in Eq. 5.6 is a Lorentz invariant. 


5.8 Formal Structure of Maxwell’s Theory 


Magnetic theory and its set of equations, in essence, now stands before us. To make 
the connection with its usual three-dimensional formulation in terms of electric field 
E and magnetic field B, we define E and B in each inertial frame by the first of the 
following equations: 
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€ —E, 0 —cB, cB, 
a —E> cB3 0 —cB, 
—E; —cB, cB 0 
so (5.123) 
O —-E, —-E, —E£; 
eu = E, 0 —cB; cB, 
E, cB, 0 —cB, 
E3; —cB. cB, 0 


which, by making the duality replacement E— cB and cB E, yields in the 
following tensor set: 


QO —cB, —cB, —cB; 


e cB, 0 —E3 E> 

oe | cBy cE3— 0 EY 
cB, —E, E, 0 

so (5.124) 

0 cB, cB> —cB3 

Ber = E, 0 —cB, cE 
E> cB3 0 —E, 
EB, —cB, cB, 0 


The two invariants of E“’—immediately recognizable as such from their mode of 
formation—can be expressed as follows: 


1 
L5G) 

1 - 3 
Y = Bo" = cB E 


eae 
1" =} a" = olf — ||? 


(5.125) 


But if we expand the first equation into the form of Stueckelberg-Lagrangian density 
in the following equation form, we have [1]: 

yec? eck? 
2 








2 
£= Fp PH + J,A" (OA) = aa) (5.126) 
Note that this equation is known as the Lagrangian density equation, or what we 
call the More Complete Electromagnetic (MCE) theory in order to generate a scalar 
longitudinal wave (SLW); more information about such matters is addressed in 
Chap. 6 of this book. 
Equation 5.126 will be written in a general form based on the following argu- 


ment. That is, the Helmholtz theorem (e.g., Griffiths, 2007 [2]) states that any three- 
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dimensional vector (electrical current density J, in this case) can be uniquely 
decomposed into circulating- and _ gradient-driven components, — thus 
j= Vet Vv x G. Here x and G are scalar and vector functions of spacetime, 
respectively. 

Using the result of Woodside [4] provides a unique decomposition for a smooth 
Minkowski field, and using only 4-vetor geometry will produce a new form of 
Maxwell’s equations as follows: 


=; OA 
E=-Vo-= (5.127) 
Ra Vea (5.128) 

- 1 0¢ 
CaVeA+ ia, (5.129) 

« LOE = & 1 0G. os 
Cc 
= » Of=s= > 10f| p 
WES LO ara seh 2 (5.131) 
Cc 


A and @ are the classical vector and scalar potentials, respectively. B and E are the 
magnetic and electric fields, respectively, as well. In the preceding set of Maxwell’s 
equations, constant C is a non-zero dynamic field with the same dimensions as the 
magnetic field. 

Bear in mind that in classical electrodynamics, C = 0 is known as the Lorentz 
gauge. Also, in the previous set of equations p and J are electric charge density and 
current charge density, respectively, while ¢ and yw represent the permittivity and 
permeability of media, not necessarily of the vacuum, respectively. 

The homogeneous equations V - B=Oand V x E+ (0B / or) = Oare identical 
to the classical model, yielding Eqs. 5.129 and 5.130. Using Eqs. 5.127 through 
5.131, the most general Lagrangian density for a mass-less, 4-vector field (A”), 
which is no more than quadratic in its variables and derivatives, is given by the 
following equation based on Eq. 2 in Stueckelberg’s 1938 article [1]—that is, 
Stueckelberg-Lagrangian density: 








2 Ie 2 2 Ot 


(5.132) 


> 2 A > 32 2 7 = i : 
eee (va +34) (Vx A) | -po4+7-4- (7-44.35) 


The result we obtained in Eq. 5.126 will be applied to generate a SLW, where we use 
it as a Lagrangian density equation for what we call the MCE theory. 
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